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1 



1 Introduction 



Riccati equation (RE) 



X 



a(x)4> 2 + b(x)4> + c(x) 



(1) 



is one of the most simple nonlinear differential equations because it is 
of first order and with quadratic nonlinearity. Obviously, this was the 
reason that as soon as Newton invented differential equations, RE was the 
first one to be investigated extensively since the end of the 17th century pQ. 
In 1726 Riccati considered the first order ODE 



with polynomial in x function u(x). Evidently, the cases degw = 1, 2 corre- 
spond to the Airy and Hermite transcendent functions, respectively. Below 
we show that Hermite transcendent is integrable in quadratures. As to Airy 
transcendent, it is only F-integrable 1 though the corresponding equation it- 
self is at the first glance a simpler one. 

Thus, new transcendents were introduced as solutions of the first order 
ODE with the quadratic nonlinearity, i.e. as solutions of REs. Some classes 
of REs are known to have general solutions, for instance: 



where all a,b,a are constant in respect to x. D. Bernoulli discovered (1724- 
25) that this RE is integrable in elementary functions if a = —2 or a = 
—Ak(2k — 1), k — 1, 2, 3, Below some general results about RE are pre- 
sented which make it widely usable for numerous applications in different 
branches of physics and mathematics. 

2 General solution of RE 

In order to show how to solve in general form, let us regard two cases. 



w x = w 2 + u(x) 



y + ay 2 = bx 



.a 



2.1 a(x)=0 

In case a(x) =0, RE takes particular form 



4> x = b(x)4> + c(x), 



(2) 



iSee Ex.3 



2 



i.e. it is a first-order LODE and its general solution can be expressed in 
quadratures. As a first step, one has to find a solution z(x) of its homogeneous 
part 2 , i.e. 

z(x) : z x = b(x)z. 

In order to find general solution of Eq.(J2J) let us introduce new variable 

4>(x) = (p(x)/z(x), i.e. z(x)(p(x) = 4>(x). Then 

(z(x)(f)(x)) x = b(x)z(x)(fi(x) + c(x), i.e. z{x)(j>(x) x = c(x), 
and it gives us general solution of Eq. (0) in quadratures 

<j)(x) = z(x)(f)(x) = z(x)( / dx + const). (3) 

J z{x) 

This method is called method of variation of constants and can be 

easily generalized for a system of first-order LODEs 

y = M x )y + f( x )- 

Naturally, for the system of n equations we need to know n particular so- 
lutions of the corresponding homogeneous system in order to use method of 
variation of constants. And this is exactly the bottle-neck of the procedure - 
in distinction with first-order LODEs which are all integrable in quadratures, 
already second-order LODEs are not. 



2.2 a(x) ^ 

In this case one known particular solution of a RE allows to con- 
struct its general solution. 

Indeed, suppose that <pi is a particular solution of Eq.((T}, then 

c = <pi, x - a<fil - htp\ 
and substitution <fi = y + (pi annihilates free term c yielding to an equation 

y x = ay 2 + by (4) 
with b = b + 2a<fi. After re- writing Eq.(JH) as 



Vx , b 

2 see Ex.1 



= a + - 

r y 
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and making an obvious change of variables <pi = l/y,we get a particular case 
of RE 

4>i,x + b<fii - a = 

and its general solution is written out explicitly in the previous subsection. 



Example 2.1 As an important illustrative example leading to many ap- 
plications in mathematical physics, let us regard a particular RE in a form 

y x + y 2 = x 2 + a- (5) 

For a — 1, particular solution can be taken as y = x and general solution 
obtained as above yields to 



y = x + 



e x2 



f e~ x2 dx + const ' 



i.e. in case (JHJ) is integrable in quadratures. Indefinite integral J e~ x ' dx 
though not expressed in elementary functions, plays important role in many 
areas from probability theory till quantum mechanics. 

For arbitrary a, Eq.((HJ) possess remarkable property, namely, after an 
elementary fraction-rational transformation 

y = X + -^— (6) 

y + x 



it takes form 



Vx + y 2 = + a, a = a + 2, 



i.e. form of original Eq.fjSJ) did not change while its rhs increased by 2. In 
particular, after this transformation Eq.fjSJ with a = 1 takes form 

y x + y 2 = x 2 + 3 

and since y = x is a particular solution of then y = x + l/xis& 
particular solution of the last equation. It means that for any 

a = 2k + l, jfe = 0,1,2, ... 

general solution of Eq. (0) can be found in quadratures as it was done for the 
case a = 1. 

In fact, it means that Eq.© is form-invariant under the transformations 
(jBJ). Further we are going to show that general RE possess similar property 
as well. 
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2.3 Transformation group 

Let us check that general fraction-rational change of variables 

~ = a(x)(f> + /3(x) 

9 7 (ar)0 + 5(x) 1 ] 

transforms one Riccati equation into the another one similar to Example 
2.1. Notice that (J7J) constitutes group of transformations generated by 

-, a(x)4>, <j> + p(x), 



thus only actions of generators have to be checked: 



• 0=1/0 transforms into 

4> x + c(x)4> 2 + b(x)4> + a(x) = 0, 

• = a(x)(f) transforms into 

4>x - ^y0 2 _ [K x ) + ( lo g"( :E ))x]0 - a(x)c(x) = 0, 

• = + (3{x) transforms Q into 

4> x — a{x)(j) 2 + [2fi(x)a(x) — b{x)]4> — c = 0, 

where 

c = a(x)p 2 (x) — b(x)/3(x) + c{x) + (3{x) x . 

Thus, having one solution of a some Riccati equation we can get immedi- 
ately general solutions of the whole family of REs obtained from the original 
one under the action of transformation group (|7j) . 

It is interesting to notice that for Riccati equation knowing any three 
solutions 0i, 02, 03 we can construct all other solutions using a very simple 
formula called cross-ratio: 

0-02 03 -02 

with an arbitrary constant A, where choice of A defines a solution. In order 
to verify this formula let us notice that system of equations 
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'0 = a(x)<p 2 + b(x)4> + c(x) 

(f>i = a{x)4>\ + b(x)(f>x + c(x) 

02 = a(x)4>l + H x )02 + c(x) 
^4> 3 = a(x)4>l + K x )4>3 + c{x) 

is consistent if 

~4> <f) 2 <f> l" 

01 0i i =0 

02 01 02 1 
.03 01 03 1. 

and direct calculation shows that this condition is equivalent to 

d , 0-01 03-0i v _ Q m 

dxV-02 03 - 02^ ' ^ ^ 

As it was shown, REs are not invariant under the action of (J7J) while 
(0) conserves the form of equations but not form of the coefficients. On the 
other hand, it is possible to construct new differential equations related to 
a given RE which will be invariant with respect to transformation group (jJJ) 
(see next section). 

At the end of this section we consider a very interesting example jHj 
showing connection of Eq. © with first integrals for generalization of one of 
Kovalevskii problems |HJ. 



Adler's example System of equations 



Vj,x + 2yj = syj, s = Y^ Vj, J = 2 > 



..n 



(10) 



was studied by Kovalevskii in case | n = 3 | and it was shown that there exist 
two quadratic first integrals 

^i = (Vi - 2/2)2/3, F 2 = (y 2 - 2/3)2/1 
and therefore Kovalevskii problem is integrable in quadratures. 

the use of (jHJ) gives us immediately following some first 

m - y% Vk - y% 

Vi - Vj Vk - Vj ' 



In case 
integrals 



n > 4 
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i.e. Svs.([10|) has nontrivial first integrals for arbitrary n. 

It is interesting that for this example solution of Svs. (jlO|) is easier to 
construct without using its first integrals. Indeed, each equation of this 
system is a Riccati equation if a is regarded as given, substitution yi = 

(f>i, x /2(f) i gives 

i,xx S&i, xxi (f)i,xx &{x) ~\~ Cj, S dxx/ &X 

and equation for a has form 

&XX 0i x / 1 ^ ^ 1 v 

a x 2 a — c\ a — c n 

After integration a 2 x = const(a — Ci)...(a — c„), i.e. problem is integrable in 
quadratures (more precisely, in hyper-elliptic functions). 



In fact, one more generalization of Kovalevskii problem can be treated 
along the same lines - case when function s is not sum of yj but some 
arbitrary function s = s(x\, ...,x n ). Then equation on a takes form 

&XX / a x &X \ 

= a x s[ h ... H 



d x a — C\ a — c n 

which concludes Adler's example. 



2.4 Singularities of solutions 

All the properties of Riccati equations which have been studied till now, are 
in the frame of local theory of differential equations. We just ignored possible 
existence of singularities of solutions regarding all its properties locally, in 
a neighborhood of a point. On the other hand, in order to study analytical 
properties of solutions, one needs to know character of singularities, behavior 
of solutions at infinity, etc. 

One can distinguish between two main types of singularities - singularities, 
not depending on initial conditions (they are called fixed) and depending on 
initial conditions (they are called movable). Simplest possible singularity 
is a pole, and that was the reason why first attempt of classification of the 
ordinary nonlinear differential equations of the first and second order, sug- 
gested by Painleve, used this type of singularities as criterium. Namely, list 
of all equations was written out, having only poles as movable singulari- 
ties (see example of PI in Chapter 1), and nice analytic properties of their 
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solutions have been found. It turned out that, in particular, Painleve equa- 
tions describe self-similar solutions of solitonic equations (i.e. equations in 
partial derivatives): P2 corresponds to KdV (Korteweg-de Vries equation), 
P4 corresponds to NLS (nonlinear Schrodinger equation) and so on. 

Using cross-ratio formula (JBJ), it is easy to demonstrate for a Riccati equa- 
tion that all singularities of the solution 0, with an exception of singu- 
larities of particular solutions <fii, $2, $3, are movable poles described as 
following: 

03 = 1 _ A (02 - Afa) 

where A is a parameter defining the solution 0. Let us construct a solution 
with poles for Eq.fjHJ) from Example 2.1. We take a solution in a form 

y— — \- a + a t (x + e) + a 2 (x + e) 2 + ... (11) 

x + e 

with indefinite coefficients a,, substitute it into (0) and make equal terms 
corresponding to the same power of (x + e). The final system of equations 
takes form 

a = 0, 

3<2i — a — e 2 = 0, 
4a 2 + 2s = 0, 
< 5a 3 — 1 + a\ = 0, 
6a 4 + 2aia 2 = 0, 
7a5 + 2aia 3 + a\ = 

1. 

and in particular for a = 3, e = the coefficients are 

a! = 1, 02 = 03 = ... = 
which corresponds to the solution 

1 

y = x H — 
x 

which was found already in Example 2.1. 

This way we have also learned that each pole of solutions have order 1. 
In general case, it is possible to prove that series (jllj) converges for arbitrary 
pair (e, a) using the connection of RE with the theory of linear equations 



8 



(see next section). In particular for fixed complex a, it means that for any 
point xq = —e there exist the only solution of flU) with a pole in this 
point. 

As to nonlinear first order differential equations (with non-quadratic non- 
linearity), they have more complicated singularities. For instance, in a simple 
example 

Vx = y 3 + 1 

if looking for a solution of the form y = ax k + .... one gets immediately 

akx*' 1 + ... = a 3 x 3k + ... k - 1 = 3k => 2/c = — 1 

which implies that singularity here is a branch point, not a pole (also see j2]). 
It make RE also very important while studying degenerations of Painleve 
transcendents. For instance, (jHJ) describes particular solutions of P4 (for 
more details see Appendix). 



3 Differential equations related to RE 
3.1 Linear equations of second order 

One of the most spectacular properties of RE is that its theory is in fact 
equivalent to the theory of second order homogeneous LODEs 

ipxx = b(x)if) x + c(x)i) (12) 

because it can easily be shown that these equations can be transformed into 
Riccati form and viceversa. Of course, this statement is only valid if Eq. (|T|) 
has non-zero coefficient a(x), a(x) ^ 0. 

► Indeed, let us regard second-order homogenous LODE (fT2"|) and make 
change of variables 



which implies 



i>x , ^xx V x 

— , then (j) x = — -, 



^ xx X i i , i2 

— = <i>x + ^ = <i>x + <i> 



and after substituting the results above into initial LODE, it takes form 

<fix = 4> 2 + b(x)(j) + c(x). 
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which is particular case of RE. I 

•4 On the other hand, let us regard general RE 

<p x = a(x)4> 2 + b(x)(f> + c(x) 

and suppose that a(x) is not = while condition of a(x) = transforms RE 
into first order linear ODE which can be solved in quadratures analogously 
to Thomas equation (see Chapter 1). Now, following change of variables 

= tl 
a(x)ib 

transforms RE into 



^xx 1 (^x\ 2 . a(x) x if) x / ip x \ 2 b(x)ifi x , . 
H — -j-r ( — J + / \n ~i~ ~~ a(x) ) --— + c(x) 



a(x)ifj a(x) \ ib ) a(x) 2 ib \a(x)ibJ a(x) ib 

and it can finally be reduced to 

a(x)ib xx — \_a(x) x + a(x)b(x)]ipx + c(x)a(x) 2 ip = 
which is second order homogeneous LODE. I 

Now, analog of the result of Section 2.2 for second order equations can 
be proved. 

Proposition 3.1 Using one solution of a second order homogeneous LODE, 
we can construct general solution as well. 

► First of all, let us prove that without loss of generality we can put 
b(x) = in if) xx + b(x)ip x + c(x)ib = 0. Indeed, change of variables 

^(x) = e ~y b( ^ dx ^(x) ib x = - -^)e~^ / b(x)dx 



and finally 



i>xx + = 0, c = c - ^b 2 - hj x . (13) 



Now, if we know one particular solution ip\ of Eq. (jlH|1 . then it follows 
from the considerations above that RE 

(j) x + (b 2 + c(x) = 
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has a solution (pi = tpi^/^i- The change of variables = (f) — 0i annihilates 
the coefficient c(x): 

(0 + 0i)' + (0 + 0i) 2 ) +c(z) = 4 + 2 + 20!0 = 



(i) a , = l + 2^ L i, (14) 
(f) cj) 

i.e. we reduced our RE to the particular case Eq.(J2J) which is integrable 
in quadratures. Particular solution z — 1/0 of homogeneous part of Eq. (jl4j) 
can be found from 

o /2 

z x = 2z— — as z —ipi 

fa 



and Eq.(j3J) yields to 



''^•'■) / ( 15 ) 



Obviously, two solutions and ?/>2 are linearly independent since Wron- 
skian <ip\,ip2 > is non- vanishing 3 : 



< i>i,ip2 >■- 



ih ^2 

$1 $2 



^1^2 ~ ^2^1 = 1^0. 



Thus their linear combination gives general solution of Eq. (|12j) .l 

Proposition 3.2 Wronskian < i/)i,ip2 > is constant iff ipi and ip 2 are so- 
lutions of 

ipxx = c(x)ip. (16) 
► Indeed, if ipi and ip2 are solutions, then 

(V>iV>2 - i>2if>'i)' = i>ii>2 - 4>24>i = c(x)(v>i^2 - V'l^) = =>- 

=>- V'i^ — ^2^1 — const . 

I 

•4 if Wronskian of two functions ip\ and ip 2 is a constant, 
ipiip 2 ~ ^2^1 — const =>- ipiip 2 ~ = 

3 see Ex.2 
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I 

Conservation of the Wronskian is one of the most important char- 
acteristics of second order differential equations and will be used further for 
construction of modified Schwarzian equation. 

To illustrate procedure described in Proposition 3.1, let us take Hermite 
equation 

XX - 2x0 x + 2X0 = 0. (17) 
Change of variables z = x /0 yields to 

_ 0xx_ _ ^ _|_ z 2 _ 2 XZ _|_ 2\ — 



and with y = z — x we get finally 

y x + y 2 = x 2 -2A-1, 

i.e. we got the equation studied in Example 2.1 with a = — 2A — 1. It means 
that all solutions of Hermite equation with positive integer A, X = n, n £ N 
can easily be found while for negative integer A one needs change of variables 
inverse to @: 



y = —x + 
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y -x 

It gives us Hermite polynomials 



(y — x)(y + x) = 7, 7 = a — 1, a — a + 2. 



'x = 


0, 


y 


X = 


1. 


y 


X = 


2. 


y 



-X, 0=1 

-x + h, 



A = n, y 



-x + 



2x 



HJx) 



as solutions. 

Notice that the same change of variables 

= — 



18) 



which linearized original RE, was also used for linearization of Thomas 
equation and Burgers equation in Chapter 1. This change of variables is 
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called log-derivative of function ip or D x \og(ip) and plays important role 
in many different aspects of integrability theory, for instance, when solving 
factorization problem. 



Theorem 3.3 Linear ordinary differential operator L of order n could be 
factorized with factor of first order , i.e. L = M o (d x — a) for some operator 
M, iff 

a = — , where ip G Ker(L). (19) 



► L = M o (d x — a), a = ip x /ip implies (d x — a)if> = 0, i.e. ip e Ker(L). I 

< Suppose that ipi — 1 is an element of the Ker(L), i.e. ipi G Ker(L). It 
leads to a = and operator L has zero free term and is therefore divisible 
by d x . 

If constant function tpi = 1 does not belong to the kernel of initial oper- 
ator, following change of variables 

ifj = — 

lead us to a new operator 

L = f- l Lof (20) 
which has a constant as a particular solution ipi for f = ip 1 . I 



Remark. Operators L and L given by (j20|) . are called equivalent oper- 
ators and their properties will be studied in detailed in the next Chapter. 

Notice that Theorem 3.3 is analogous to the Bezout's theorem on divis- 
ibility criterium of a polynomial: A polynomial P(z) =0 is divisible on the 
linear factor, P(z) = P\(z)(z — a), iff a is a root of a given polynomial, i.e. 
P{a) = 0. Thus, in fact this classical theorem constructs one to one corre- 
spondence between factorizability and solvability of L{tj)) = 0. 

The factorization of differential operators is in itself a very interesting 
problem which we are going to discuss in details in Chapter 3. Here we will 
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only regard one very simple example - LODO with constant coefficients 

L(V0 :=_ + a 1 ^ =r + ... + ^ = 0. 
In this case each root Aj of a characteristic polynomial 

A n + aiA n_1 + . . . + a n = 
generates a corresponding first order factor with 

\ _ Vfc 

Aj — — 

ant it yields to 

ip x = Xiijj => ip = de XlX 

and finally 

-_. d di / d \ \ / d> -x \ 

L = — + ai———^ + . . . + a n = (— - Ai) • • • (— - A„). 
ax n ax n 1 ax ax 

This formula allows us to construct general solution for L(ip) = 0, i.e. for 
ip G Ker(L), of the form 

^ = J2 CieXiX 

in the case of all distinct roots of characteristic polynomial. 

In case of double roots A& with multiplicity it can be shown that 

^ = J2Pk(x)e x « x (21) 

where degree of a polynomial Pk(x) depends on the multiplicity of a root, 
degPk(x) < m k — 1 (cf. Ex.4) 

3.2 Schwarzian equation 

Let us regard again second-order LODE 

V>x* + 6(^)V>x + c(x)i) = (22) 

and suppose we have two solutions 4>i,ip2 of (}2"2"j) . Let us introduce new 
function </? = ip 1 /ip 2 , then 
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V ' = J* ' ^ = 6 ^ V| + 2 ^ 



which yields 



Vxx u , that 

= -b{x) - 2— 



and substituting = ^ = (log ^2)2 into (pTJ) related to ((221) we S e ^ finally 

' fl «y_i^ = ( x (23) 

Left hand of (|23|) is called Schwarz derivative or just Schwarzian and 

is invariant in respect to transformation group (|IJ) with constant coefficients 

acp + (3 



7V? + 5 

It is sufficient to check only two cases: 



ip = — and = aip + (3. 



which can be done directly. 



This equation plays major role in the theory of conform transformations 
of polygons [Zj . 



3.3 Modified Schwarzian equation 

Notice that substitution (p = ipi/ip2 allowed us to get invariant form of the 
initial Eq. (j22j) . Another substitution, namely, tp = ip\ip2i leads to similar 
equation which differs from classical Schwarzian equation (|23|) only by a 
constant term and this is the reason why we call it modified Schwarzian 
equation. This form of Schwarzian equation turns out to be useful for a 
construction of approximate solutions of Riccati equations with parameter 
(see next section). In order to construct modified Schwarzian equation, we 
need following Lemma. 

Lemma 3.4 Let ipi, ^2 are two linear independent solutions of 

ipxx = c(x)ip. (24) 

Then functions 
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constitute a basis in the solution space of the following third order equation: 

Vxxx = 4c(x)ip x + 2c x (x)<p. (25) 

► Using notations 

<Pi = iPi, V2 = ipl, <pi = 1P11P2, 
we can compute Wronskian W of these three functions 

W =< (fl, <f 2 , ¥3 >= {^>l^2,x ~ 1p2lpl,x) 3 =< 01, 02 > 3 

and use Proposition 3.2 to demonstrate that 

W = const ^ 0, 

i.e. functions (pi are linearly independent. 
After introducing notations 

V=< -01,-02 > and ./• 



it is easy to obtain 



— = /2"/l, ^ = / 2 + /l 

V?3 V?3 



which yields to 

Substitution of these fj into 



fl = ^y., a = *^. (26) 

2<^3 2v?3 



A> + /J = c(x) 

gives 

4c(a;)^ 2 + <pl - 2ipip xx = V 2 . (27) 

with ip = ip 3 and differentiation of Eq.(j2ZJ) with respect to x gives Eq. (|23J) 
and it easy to see that equations (J27|) and (1251) are equivalent. 

Analogous reasoning shows that ipi, tp 2 are also solutions of Eq.(|23|). I 

Equation ({25} as well as its equivalent form (|2*7|) will be used further for 
construction of approximate solutions of REs, they also define solitonic hier- 
archies for KdV and NLS. It will be more convenient to use (|2*7jl in slightly 
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different form. 



Let us rewrite (l2?1) as 



4cW + 4_!^ 

Lf z if Lf z 



and introduce notation a = 1/ip, then 

3 a? la 



c(x) = r 3 _ tT + vV (» 



and compare this equation with Schwarzian equation 

/ \ 2 
/ \ 3 / lf xx \ I (p x 

c[x) — - 



4 \(p x ) 2 ip a 



one can see immediately why Eq. (j28|) is called modified Schwarzian 
equation. 

Notice that after the substitution a = e 2b , rhs of modified Schwarzian 
equation, i.e. modified Schwarzian derivative, Dmod, takes a very 
simple form 

Dmod(a) :=- A ^ 2 - l - a — = b xx + bl 
4 a z 2 a 

which is in a sense similar to -D x log. Indeed, for ifj = e 9 , 

D x log(V>) = ^ = if x = e ~*-^fi 

while 

d 2 

Dmod(e 2lp ) = e-*—e v . 

dx 2 



At the end of this section let us stress the following basic fact: we have 
shown that from some very logical point of view, first order nonlinear Riccati 
equation, second order linear equation and third order nonlinear Schwarzian 
equation are equivalent! It gives us freedom to choose the form of equation 
which is most adequate for specific problem to be solved. 
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4 Asymptotic solutions 



In our previous sections we have studied Riccati equation and its modifica- 
tions as classical ordinary differential equations, with one independent vari- 
able. But many important applications of second order differential equations 
consist some additional parameter A, for instance one of the most significant 
equations of one-dimensional quantum mechanics takes one of two forms 

= (X + u)ij (29) 
^xx = (X 2 + u 1 X + u 2 )^ (30) 

where Eq. ()2T)|) is called Schrodinger equation and Eq. (f3T)|) can be consid- 
ered as modified Dirac equation in quantum mechanics while in applications 
to solitonic hierarchies it is called Zakharov-Shabat equation. Notice that 
Schrodinger equation with u(x) = x 2 equivalent to Eq. (fTT|) . Coefficients of 
these two equations have special names - u,Ui, 112 are called potentials due 
to many physical applications and A is called spectral parameter because 
of following reason. Schrodinger equation, being rewritten as 

L(ip) = \ip, L(ip) = %p xx - uip 

becomes obviously an equation for eigenfunctions of operator L (with appro- 
priate boundary conditions, of course). This operator is called Schrodinger 
operator. 

For our convenience we name the whole coefficient before ip as general- 
ized potential allowing it sometimes to be a polynomial in A of any finite 
degree. Coming back to Eq.(J2U), the generalized potential is just the func- 
tion c(x). 

Now, with the equation having a parameter, problem of its integrability 
became, of course, more complicated and different approaches can be used 
to solve it. If we are interested in a solution for all possible values of a 
parameter A, asymptotic solution presented by a formal series can always 
be obtained (section 4.1) while for some specific exact solutions can be con- 
structed (sections 4.2, 4.3) in a case of truncated series. It becomes possible 
while existence of a parameter gives us one more degree of freedom to play 
with. Cf. with Example 4.2 where exact solution has been obtained also as 
a series and its convergence resulted from the main theorem of the theory of 
differential equations on solvability of Cauchy problem. On the other hand, 
this solution is valid only for some restricted set of parameter ' values, namely 
for integer odd a. 
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4.1 RE with a parameter A 

Let us show first that the RE with a parameter A corresponding to Eq. (}30|) . 
namely 

U + f = X 2 + u 1 X + u 2 , with f = D x \og(^), (31) 
has a solution being represented as a formal series. 



Lemma 4.1 Eq.([31j) has a solution 



/ = A + / + y + -.. 
where coefficients /,■ are differential polynomials in u\ and u 2 - 



(32) 



► After direct substituting the series ({HZ} into the equation for / and 
making equal corresponding coefficients in front of the same powers of A, we 
get 



2/o = «i 

2/l + / ,x + / 2 = «2 

2/ 2 + + 2/0A = 

2/3 + / 2 ,x + 2/o/ 2 + ft = 



and therefore, coefficients of 
and u%. I 



are differential polynomials of potentials U\ 



Notice that taking a series 
9 = 



9i , 92 , 93 



X + 90 + J + V + V 



(33) 



as a form of solution , we will get a different system of equations on its co- 
efficients gf. 



-2g = ui 

-2g t + g 0>x + gl = u 2 
-2# 2 + 9i, x + 2g gi = 
-2c/ 3 + 92,x + 2go92 + 9i=0 
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Solution of Eq.([3ip constructed in Lemma 4.1 yields to the solution of original 
Zakharov-Shabat equation ()30|) of the form 

A) = = e^( Vo (x) + + + ^ + ....) (34) 

and analogously, the second solution is 

M*, A) = J*"** = e-^(Ux) + ^ + ^# + ^# + ..••) (35) 

A A z A 13 

In fact, it can be proven that Wronskian < ipi, -0 2 > is a power series on A (see 
Ex.8) with constant coefficients. Notice that existence of these two solutions 
is not enough to construct general solution of initial Eq.(j30J) because linear 
combination of these formal series is not defined, also convergence problem 
has to be considered. On the other hand, existence of Wronskian in a conve- 
nient form allows us to construct family of potentials giving convergent series 
for ()34|) and (J33j) . We demonstrate it at the more simple example, namely 
Schrodinger equation (j2HJ). 

Let us regard Schrodinger equation ()29j) . its RE has form 

f x + f = X + u, with / = D a! log(V), (36) 

and it can be regarded as particular case of (JHDJ), i.e. the series for / yields 
to 

f = k + f + j + ... } \ = k\ (37) 

and g(x, k) = f(x, —k). We see that in case of ()29j) there exists a simple way 
to calculate function g knowing function / and it allows us to construct two 
solutions of Schrodinger equation (J2HJ): 

Mx, k) = eSf^ = e k *(l + ^ + ^ + + - (38) 

and 

Substitution of say ipi into (}2T?j) gives a recurrent relation between coefficients 

Cj+i,x = \{ U Q - Cj,xx), Co = 1- (39) 

In particular, 

u = 2(i, x (40) 
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which means that in order to compute potential u it is enough to know 
only one coefficient C\ of the formal series! Below we demonstrate how 
this recurrent relation helps us to define potentials corresponding to a given 
solution. 



Example 4.2 Let us regard truncated series corresponding to the solutions 
of (El 

^i = e^(l + |), ^ 2 = e-^(l-|), 

then due to flHUJ) 

u — 2(i,x, (i,xx = 2Ci,xCi 
and Wronskian W of these two functions has form 

W =< Vi, ^2 >= -2k + i(Ci - Cx,x). (41) 

Notice that 

(Ci - Ci,x) x = o 

and it means that W does not depend on x, W = W(A;). Introducing notation 
ki for a zero of the Wronskian, W(/ci) = 0, it is easy to see that 

d 2 - Ci,, = kl 

which implies that ipi and ip2 are solutions of (|29j) with 



^ = h ~ l + e^Mx-xo) = taIlh ^ ~ 

and potential 

(2h) 2 -2k\ , s 

u = —2 i — = - (42) 

( e ki(x-x ) + e -k x (x-x )y cosh 2 (A;i(a;-xo))' 

where x is a constant of integration. 

It is important to understand here that general solution of Schrodinger 
equation ()29|) can be now found as a linear combination of ipi and ip2 for all 
values of a parameter A = k 2 with exception of two special cases: k = 
and k = ki which implies functions ip\ and ip2 are linearly dependent in 
these points. 

Fig.l (...) 
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At the Fig. 1 graph of potential u is shown and it is easy to see that 
magnitude of the potential in the point of extremum is defined by zeros of the 
Wronskian W. At the end of this Chapter it will be shown that this potential 
represents a solitonic solution of stationary KdV equation, i.e. solution of 
a Riccati equation generates solitons! 

4.2 Soliton-like potentials 

It this section we regard only Schrodinger equation ()29)1 and demonstrate 
that generalization of the Example 4.2 allows us to describe a very impor- 
tant special class of potentials having solutions in a form of truncated series. 

Definition 4.3 Smooth real-valued function u(x) such that 

u(x) ^ for x — > ±oo, 

is called transparent potential if there exist solutions of Schrodinger equa- 
tion ()29|) in a form of truncated series with potential u(x). 

Another name for a transparent potential is soliton-like or solitonic 
potential due to many reasons. The simplest of them is just its form which 
is a bell-like one and "wave" of this form was called a soliton by [10 and this 
notion became one of the most important in the modern nonlinear physics, 
in particular while many nonlinear equations have solitonic solutions. 

Notice that truncated series ipi and ip2 can be regarded as polynomials 
in k of some degree N multiplied by some exponent (in Example 4.2 we had 
N — 1). In particular, it means that Wronskian W =< ipi,ip2 > is odd 
function, W(— k) = —W(k), vanishing at k = and also it is a polynomial 
in k of degree 2N + 1: 

N 

W{k) = -2k\[{k 2 -k]). (43) 
i 

As in Example 4.2, functions tp\ and ip2 are linearly dependent at the 
points kj, i.e. 

ifa.(x,kj) = Aj^,{x,kj), j = 1,2, JV 
with some constant proportionality coefficients Aj. 
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Theorem 4.4 Suppose we have two sets of real positive numbers 

{%}, {Bj}, j = 1, 2, N, kj, Bj > 0, kj, Bj e E 

such that numbers kj are ordered in following way 

ki > k 2 > ... > k N > 

and Bj are arbitrary. Let functions i[>i(x, k),Tp 2 {x, k) have form 

rk(x, k) = e kx (k N + a x k N - 1 + ... + a N ), fo(x, k) = {-l) N ijj x {x, —k) 

with indefinite real coefficients aj,j = 1, 2, N. 

Then there exist unique set of numbers {aj} such that two functions 

■01 (x, fey), 02 (aJ, %) satisfy system of equations 

2 (x, = (-iy +1 Bj^(x, kj) (44) 
and Wronskian of these two functions has form 

N 

W(k) = -2k\[(k 2 -k 2 ). (45) 
i 

► Our first step is to construct aj. It is easy to see that (jUj) is equivalent 
to the following system of equations on {aj} (for simplicity the system is 
written out for a case N = 3) 

k\ai + E x k\a 2 + a 3 + kfEi = 

k\E 2 a x + k 2 a 2 + E 2 a 3 + fcf = (46) 

fcfai + E 3 k 3 a 2 + CL3 + k%E 3 = 

where following notations has been used: 

£j = = — = tanliTj, Tj = kjX + (3j, Bj = e 2f3j . 

(To show this, it is enough to write out explicitly ipi and ip 2 in roots of 
polynomial and regard two cases: N is odd and N is even. For instance, if 
k = ki and N is odd, we get 

0i (x, k) = e kx {k N + a^- 1 + ... + a N ) = ...0 2 (x, k) = {-lf^x, —k)... 

) 
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Obviously, < E 3 < 1 and for Ej = 1, Vj = 1,2,3 the system Ijflty 

takes form 

kfai + k x a 2 + a 3 + kf = 

fc|ai + fc 2 02 + «3 + k\ = (47) 
fcgOi + k 3 a 2 + a 3 + fcf = 

(...). Thus, it was shown that determinant of Svs. (j46j) is non-zero, i.e. all 
dj are uniquely defined and functions ipi(x, kj), ip 2 (x, kj) are polynomials. 

In order to compute the Wronskian W of these two functions, notice first 
that W is a polynomial with leading term — 2k 2N+1 . Condition of propor- 
tionality (|4^j) for functions ipi(x,kj), ^(x, fcj) provides that kj are zeros of 
the Wronskian and that W is an odd function on k, i.e. (I45J) is proven. I 

In order to illustrate how to use this theorem for construction of exact 
solutions with transparent potentials let us address two cases: N = 1 and 
N = 2. 

Example 4.5 In case N — 1 explicit form of functions 

fa = e kx (k + ai), ip 2 = e~ kx (k - a x ) 

allows us to find a x immediately: 

ai = —k\E\ = —ki tanhyi = —k\ tanh {k\x + (3i) 

which coincides with formula for a solution of the same equation obtained in 
Example 4.2 

(i = —k\ tanh k±(x — Xq) 

for xq = (3i/k%. As to potential u, it can be computed as above using recur- 
rent relation which keeps true for all N. 

Notice that using this approach we have found solution of Schrodinger 
equation by pure algebraic means while in Example 4.2 we had to solve 
Riccati equation in order to compute coefficients of the corresponding trun- 
cated series. 

The system f!4fij) for case N = 2 takes form 

Mi + + k\E x = 
E 2 k 2 ai + a 2 + k\ = 
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which yields to 



(k 2 2 - hfiEh 



h - k 2 E 1 E 2 

= D x log ({k 2 - fci) cosh(ri + t 2 ) + (k 2 + h) cos1i(ti - r 2 )) (48) 
and corresponding potential u = 2ai )X has explicit form 

u = 2D 2 X log {(k 2 - k x ) cosh(ri + r 2 ) + (A; 2 + fci) cosh(Yi - r 2 )) (49) 
where 

x — > ±00 ai — * ±(fei + fc 2 ), 

i.e. m is a smooth function such that 

it — > for x — > ±00. 



Formulae (|48jl and (|4T?|l have been generalized for the case of arbitrary iV 
by Hirota whose work gave a rise to a huge amount of papers dealing with 
construction of soliton-like solutions for nonlinear differential equations 
because some simple trick allows to add new variables in these formulae (see 
[E] and bibliography herein). For instance, if we take 

Tj — — kjX ~\~ j3j — — kjX ~\~ ^j^ 1 

then formula ()49|) gives particular solutions of Kadomtzev-Petviashvili (KP) 
equation 

{—Aut + u xxx + Quu x ) x + 3u yy = (50) 
which is important model equation in the theory of surface waves. 

4.3 Finite- gap potentials 

In previous section it was shown how to construct integrable cases of Schrodinger 
equation with soliton-like potentials vanishing at infinity. Obvious - but not 
at all a trivial - next step is to generalize these results for construction of 
integrable cases for Schrodinger equation with periodic potentials. In the 
pioneering work jTTj the finite-gap potentials were introduced and described 
in terms of their spectral properties but deep discussion of spectral theory 
lays beyond the scope of this book (for exhaustive review see, for instance, 
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[12J). The bottleneck of present theory of finite-gap potentials is following: 
spectral properties formulated by Novikov's school provide only almost pe- 
riodic potentials but do not guarantee periodic ones in all the cases. 

We are going to present here some simple introductory results about 
finite-gap potentials and discuss a couple of examples. For this purpose, most 
of the technique demonstrated in the previous section can be used though 
as an auxiliary equation we will use not Riccati equation but its equivalent 
form, modified Schwarzian (|28)1. 

Generalization of Lemma 4.1 for the case of arbitrary finite polynomial 
c(x) can be formulated as follows. 



Lemma 4.6 Equation for modified Schwarzian 

7lf - + Xmk2 = U< < X > A ) : = ^ + U ^ + - + «m (51) 

4 hr 2 h 

with any polynomial generalized potential U (x, A) has unique asymptotic 
solution represented by formal Laurant series such that: 



h(x,X) = 1 + J2 X ~ kh k(x) (52) 

k=\ 

where coefficients hj are differential polynomials in all U\, ...,u m . 

► The proof can be carried out directly along the same lines as for Lemma 
4.1. I 

Direct corollary of this lemma is following: coefficients of the formal so- 
lution are explicit functions of potential U(x,X). In particular, for m = 1 
which corresponds to Schrodinger equation ()29)1 with generalized polynomial 
potential A + u we have 

hi = -u, 2h 2 = -hi, xx - hi, ... (53) 

Definition 4.7 Generalized potential 

U(x,X) = X m + u 1 X m - l + ... + u m 

of an equation 

?Pxx = U(x, X)ip 
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is called iV-phase potential if Eq. (J2oTl . 

^xxx 

= 4U(x,X)<p x + 2U x (x)ip, 
has a solution (p which is a polynomial in A of degree N: 

N 

<p( x , X) = X N + (p l (x)X N ~ 1 + ... + VN (x) = J](A - 7j (x)). (54) 

3=1 

Roots jj(x) of the solution (p(x, A) are called root variables. 

In particular case of Schrodinger equation this potential is also called 
finite-gap potential. As it follows from ^2], original "spectral" definition 
of a finite- gap potential is equivalent to our Def. 4.7 which is more conve- 
nient due to its applicability not only for Schrodinger equation but also for 
arbitrary equation of the second order. 

There exists direct connection between the statement of Lemma 4.6 and 
the notion of finite-gap potential. One can check directly that if function 
h(x, A) is solution of Eq. lfBTj) . then function <p(x, A) = l/h(x, A) is solution of 
Eq. (j25|) and can be written as formal series 

oo 

(p(x,X) = l + ^2\- k <p k (x) (55) 

k=l 

with coefficients which are explicit functions of generalized potential U(x, A). 
In case when series f!55|) becomes a finite sum, we get finite-gap potential. 

Example 4.8 Let a solution tp(x, A) = A — 7 (x) is a polynomial of first 
degree and potential is also linear, i.e. m = 1, N = 1. Then after integrating 
the equation from definition above, we get 

4(A + u)</ + <f 2 x - 2ipip xx = c(A) (56) 

with some constant of integration c(A) and left part of ()56j) is a polynomial 
in A of degree 3, 

C(X) = 4A 3 + Cl A 2 + c 2 A + c 3 = 4(A - Ai)(A - A 2 )(A - A 3 ), (57) 

where Ai are all roots of the polynomial C(A). Eq. ()56j) is identity on A and 
therefore without loss of generality we write further C(A) for both sides of 
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it. This identity has to keep true for all values of A, in particular, also for 
A = j(x) which gives 

7 2 = C( T ) = 4( 7 - A0( 7 - A 2 )( 7 - A 3 ). (58) 

Now instead of solving Eq.(|55)l. we have to solve Eq. (J55Jl which is inte- 
grable in elliptic functions. 

If we are interested in real solutions without singularities, we have to 
think about initial data for Eq.([5S]). For instance, supposing that all Xj are 
real, without loss of generality 

Ai > A2 > A3, 

and for initial data (x ,j ) satisfying 

V(x ,7o) : A 3 < 70 < A 2 , 

Eq.()58p has real smooth periodic solution expressed in elliptic functions 

u = 2 7 - X 1 - A 2 - A 3 (59) 

It is our finite-gap potential (1-phase potential) and its period can be com- 
puted explicitly as 

T _ r x ' 2 dX 

~L V / (A-A 1 )(A-A 2 )(A-A 3 )' 

We have regarded in Example 4.8 particular case mN = 1. Notice that 
in general case of mN > 1 following the same reasoning, after integration we 
get polynomial C(A) of degree 2N + m 

4U(x, AV 2 + (fl - 2(pcp xx = (7(A) := AX 2N+m + ... (60) 

and correspondingly a system of 2N + m — 1 equations on iV functions, i.e. 
the system will be over-determined. On the other hand, choice of A = jj 
makes it possible to get a closed subsystem of N equations for iV functions 
as above: 

7^ = C( 7i )/n(7;-7*) 2 - (61) 

Following Lemma shows that this over-determined system of equations has 
unique solution which is defined by Svs.(ffiT|). 
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Dubrovin's Lemma Let system of differential equations (JoTj) on root 
variables 7^ defined by (p>^|) is given with 



C(A) = 4A 2JV+m + ..., 

then following keeps true: 

f. C(A)| a=7j =^(x,A)| a=7j , j = l,...,N, 

2. expression 

ip {2(p xx H J 

is a polynomial in A of degree m and leading coefficient 1. 

M Let us notice that 

JV N-l 

A =7j 



( n^ A _ ^k)) x \\= l3 = ( - 71,0= n^ A ~ ^ n (a - 7;)) 

j=2 j=l 



= -7i,«c n(Ti - 7k) 

which implies 



^|a= 7 = -7j> JJ(7j - 7*) => C(A)| A =^ = <pl(x, A)| a=7j , 

3^k 



i.e. first statement of the lemma is proven. 



Below it will be shown that also transparent potentials themselves 

can be computed algebraically. (...) 



5 Summary 

In this Chapter, using Riccati equation as our main example, we tried to 
demonstrate at least some of the ideas and notions introduced in Chapter 
1 - integrability in quadratures, conservation laws, etc. Regarding trans- 
formation group and singularities of solutions for RE, we constructed some 
equivalent forms of Riccati equation. We also compared three different ap- 
proaches to the solutions of Riccati equation and its equivalent forms. The 
classical form of RE allowed us to construct easily asymptotic solutions repre- 
sented by formal series. Linear equation of the second order turned out to be 
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more convenient to describe finite-gap potentials for exact solitonic solutions 
which would be a much more complicated task for a RE itself while general- 
ization of soliton-like potentials to finite-gap potentials demanded modified 
Schwarzian equation. 

In our next Chapter we will show that modified Schwarzian equation also 
plays important role in the construction of a differential operator commuting 
with a given one while existence of commuting operators allows us to obtain 
examples of hierarchies for solitonic equations using Lemma 4.6. In partic- 
ular, for m — 1 coefficients hk{x) of Eq. (|52jl describe a set of conservation 
laws for Korteweg-de Vries equation (KdV) 



2. Deduce formula (|T3|) regarding < ipi,ip 2 >— 1 as a linear first order 
equation on ip 2 . 

3. Prove that for L = its kernel consists of all polynomials of degree 
< m — 1. 

4. Let functions Ax and A 2 are two solutions of ([25)1 . Prove that the Wron- 
skian < A±, A 2 >= A\A 2 — A 2 A[ is solution as well. 

5. Let the function A satisfies (12 7|) . Prove that the functions 



u t + 6uu x + u. 



XXX 



0. 



(62) 



6 Exercises for Chapter 2 

1. Prove that general solution of z = a(x)z has a form 

z(x) = e f a{x)dx . 



f ± = -D\ogA± 



A 



satisfies the Riccati equations (JH). 



6. Proof that 




2 a 



1 a. 



XX 



k 2 & a = (e ie kx + e 2 e~ kx )- 2 
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